Spin liquid correlations in Nd-langasite anisotropic Kagome antiferromagnet 
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Dynamical magnetic correlations in the geometrically frustrated NdsGasSiOw compound were 
probed by inelastic neutron scattering on a single crystal. A scattering signal with a ring shape 
distribution in reciprocal space and unprecedented dispersive features was discovered. Comparison 
with calculated static magnetic scattering from models of correlated spins suggests that the observed 
phase is a spin liquid inherent to an antiferromagnetic kagome-like lattice of anisotropic Nd moments. 

PACS numbers: 75.40.Gb,75.50.Lk,75.10.Hk 
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Geometrically frustrated magnets are at the forefront 
of research because of their diverse exotic ground states 
but also because the concept of frustration plays a ma- 
jor role in fields beyond magnetism One prototype is 
the kagome antiferromagnet, a two dimensional lattice of 
corner sharing triangles of antiferromagnetically interact- 
ing spins, which attracts strong attention as a promising 
candidate for possessing a disordered ground state. At 
the quantum limit of spins 1/2, there are converging ar- 
guments towards the stabilization of a type II spin liquid 
2] consisting in a disordered ground state breaking nei- 
ther spin nor lattice symmetry, and realizing a quantum 
spin liquid state. Its search traces back to the suggestion 
of a resonating valence bond state in the triangular lat- 
tice Heisenberg model |3j ■ At the classical limit of large 
spins the basic signature of the geometric frustration is 
the infinite degeneracy of the ground state. A disorder 
phase with short range magnetic correlations, transpos- 
able to the superposition of ground state spin configura- 
tions, can be expected, with an intrinsic dynamical na- 
ture associated with zero-point motions, localized or ex- 
tended 0,0, IE 00- Subtle differences exist according to 
the spin degrees of freedom. With planar XY or isotropic 
Heisenberg spins, because of the continuous extensive en- 
tropy, thermal fluctuations play a crucial role in selecting 
a subspace within the ground state manifold [|| through 
the order out of disorder mechanism Whether non 
linear effects may even push forward this selection to a 
unique ground state is still an open question. With Ising 
spins, a cooperative paramagnet is predicted at all tem- 
peratures with a discrete extensive ground state degen- 
eracy and exponentially decaying spin-spin correlations 
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A disorder phase is experimentally easily destabi- 
lized by any kind of small perturbation Q], like sec- 
ond order interactions, structural imperfections or off- 
stoichiometry. Accordingly, most of the few materializa- 
tions of kagome antiferromagnets available to date do or- 
der at low enough temperatures, though often into uncon- 



ventional magnetic states (non collincar, non coplanar, 
multi-propagating or incommensurate orderings, disor- 
der free spin glassiness, ...) that are in some cases sus- 
taining large magnetic fluctuations. Of interest is that 
these compounds (for instance the jarosites ^lj, the Cr 
kagome bilayers 01 > the volborthite 0] and the kagome- 
staircase compounds |l4|) are best described by Heisen- 
berg spin models and that up to now, there was no ex- 
ample of a kagome antiferromagnet with a strong magne- 
tocrystalline anisotropy calling for an Ising or a XY de- 
scription. Another point was the absence of large enough 
single crystals for three-axis inelastic neutron scattering 
in all cases of kagome antiferromagnets in which a spin 
liquid state was suspected. This unique technique allows 
one to probe the dynamical magnetic correlations both 
in space and time and was earlier successfully used in ev- 
idencing signatures of spin liquid phases in compounds 
containing a 3D pyrochlore net of m agne tic moments : 
the itinerant-electron Yo.97Sco.o3Mn2 and the insu- 
lating ZnCr 2 4 0. 

Nd-langasite, NdaGasSiOu, crystallizes in the trigonal 
space group P321 and belongs to a family of materials 
that attracted strong interest for showing better piezo- 
electric properties than quartz or lithium niobate and 
tantalate njj- We recently discovered that some mem- 
bers of this family could be also of interest as geometri- 
cally frustrated magnets |18j ■ The Nd 3+ magnetic ions in 
Nd3Ga5SiOi4 are organized in corner sharing triangles in 
well separated planes perpendicular to the c-axis. Within 
these, the Nd 3+ form a distorted kagome lattice, which 
is topologically equivalent to the ideal kagome when only 
the shortest atom bridging interactions are considered. 
We emphasize that this Nd net is fully occupied and that 
large size single crystals can be grown. 

Powders of Nd-langasite were prepared by solid state 
reactions of stoichiometric amounts of high purity oxides 
at 1420°C in air whereas single crystals were grown from 
sintered powders by the floating zone method using a mir- 
ror furnace from "Cyberstar" company. Single-crystal 
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magnetization measurements confirmed that a geometric 
frustration does materialize in this compound and re- 
vealed a strong magnetocrystalline anisotropy associated 
with crystal field effects on the ground multiplet of the 
4f3 ions (total angular momentum J = 9/2) [lq. A 

quantitative analysis of the magnetic susceptibility x a t 
high temperatures T yielded a Curie- Weiss temperature 
9 = -52 K accounting for antiferromagnetic interactions 
between the Nd moments and quadrupolar crystal field 
parameters indicating that these moments are coplanar 
rotators in the kagome planes . A change in the mag- 
netocrystalline anisotropy occurs around 33 K due to the 
higher order terms of the crystal electric field, the c-axis, 
perpendicular to the kagome planes, becoming the mag- 
netization axis at low temperature |18| . No phase tran- 
sition towards a long-range magnetic order was detected 
down to 1.6 K, a temperature much smaller than the ex- 
change coupling energy scale given by 6, neither in these 
magnetization measurements nor by elastic neutron scat- 
tering on powders, clearly suggesting that a sp in liquid 
phase could be stabilized in the compound |18j . 

Inelastic neutron scattering measurements were per- 
formed at the high flux reactor of the Institut Laue 
Langevin to probe the low temperature dynamical mag- 
netic correlation. Experiments with powders on the IN5 
time of flight spectrometer informed essentially about the 
low levels of the energy spectrum of the Nd 3+ ions as split 
by the crystal electric field. An additional much weaker 
signal was also detected around 1 meV [19J . 

In this letter, we report on the results from an inelas- 
tic neutron scattering experiment on a large single crys- 
tal, 40 mm long and 5 mm in diameter, performed on 
the cold neutron three-axis spectrometer IN14 to probe 
the weak magnetic signal around 1 meV and its distri- 
bution in the reciprocal space. The single crystal was 
mounted in a cryostat with the [100] and [010] axes in 
the horizontal scattering plane. Its temperature was 
maintained at 2 K. The neutron wavelength was se- 
lected from a vertically curved graphite monochromator 
PG(002). We used a fixed final energy of 4.66 meV with 
energy resolution around 165 /ieV. In the following, the 
wave- vectors Q are expressed via their coordinates (£,C,0) 
in reciprocal lattice units (r.l.u.) or via their modulus 
\Q\/\a*\ = v/(£ + l/2C) 2 + 3/4C 2 with |o*| = 0.899 A" 1 . 

Scans in neutron energy transfer (E-scans) at different 
positions along the [100] and [210] directions are shown 
in fig. [J At Q = (1.5 0) (along [100]), a small signal 
is visible around 0.9 meV. At Q = (1.3 -0.65 0) (along 
[210]) a small signal is also present at a slightly smaller 
energy and mixed with the tail of the incoherent peak. 
Along each direction, the intensity of the signal decreases 
at larger Q values, as expected for magnetic excitations, 
which are multiplied by the square of the magnetic form 
factor 21]. Phonon excitations would have shown a rise 
of intensity with Q 2 . A half width at half maximum 
(HWHM) T = 0.47 meV is deduced from the fit of the 
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FIG. 1: Top : Intersection of the Brillouin zones with the 
(a*, b*) scattering plane (dark gray lines). The white arrows 
indicate the Q-scans performed in both the [100] and [210] 
directions. Bottom : E-scans at 2 K for two different Q po- 
sitions in each direction (crosses in the top panel). The thin 
lines are Gaussian fits of the incoherent peak. The thick line 
is a lorentzian fit (see text). 



energy response at Q=(1.5 0), with a lorentzian times 
the detailed balance factor, which corresponds to a life 
time r = 1.4 10~ 12 s for the spin correlations. 

Scans in neutron wave vector transfer (Q-scans), along 
the [100] and [210 Jdirections, at the 0.5, 0.675, 0.85 and 
1 meV energy transfers, show a weak, gradual and global 
evolution (cf. fig. yj. The signal is minimum around 
IQI/I a *l = 2, between a broad maximum and (in most 
scans) a second weaker maximum. The first maximum is 
also observed in the intermediate direction [410] at the 
same Q value than in the other two directions, confirm- 
ing the peculiar ring shape of this signal [l9j. Its am- 
plitude decreases and its position shifts from |Q|/|a*| = 
1 to 1.5 with increasing energy. This indicates a signifi- 
cant dispersion with Q of the magnetic scattering in the 
0.5-1 meV range. The position and amplitude of a sec- 
ond maximum, when present, also vary with energy. The 
lines in fig. 2 show the results of fits to two lorentzians of 
same width, multiplied by the square of the Nd 3+ mag- 
netic form factor j^- The fitted HWHM k ps 0.49 A" 1 
corresponds to a correlation length A ~ 2 A, smaller than 
the distance between the nearest neighbors Nd 3+ ions 
(4.2 A). 
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FIG. 2: Q-scans at 2 K at several energy transfers, along the 
[100] and [210] directions. The dotted lines indicate the posi- 
tion of the first maximum. The solid lines following the mea- 
sured Q distributions (circles) are lorentzian fits (see text). 
At bottom, static scattering, calculated using first neighbors 
magnetic correlations only (solid line) and in a model involv- 
ing further neighbors correlations (dashed line), are plotted 
for the same directions for comparison (see text). 



An interpretation of these results can be proposed from 
the comparison of the measured magnetic signal with cal- 
culated neutron static scattering intensities using differ- 
ent models of antiferromagnetically interacting classical 




FIG. 3: Calculated static scattering map for the antiferromag- 
netic kagome lattice (left) and for the distorted Nd-langasite 
kagome lattice (right) using first neighbor spin correlations 
(see text). The brightness increases with intensity. 



spins on the distorted kagome lattice formed by the Nd 3+ 
ions. In all cases, the calculated Fourier-transformed spin 
correlation functions have been multiplied by the square 
of the Nd 3+ magnetic form factor 20] . The simplest 
model is obtained considering correlations between first 
neighboring spins only and spin configurations minimiz- 
ing the energy on each triangle of the structure. These 
are given as having the sum of three spins equals to zero 
for Heisenberg or XY spins and two spins in one state and 
the third spin in the opposite state for the Ising spins. 
Surprisingly, the resulting Q distribution of magnetic in- 
tensity does not depend on the spin degrees of freedom. 
This validates the interpretation of the experiment with 
this model albeit the low temperature anisotropy of the 
compound is not quantitatively determined at present. It 
is shown in fig. [21 where it is compared to the same cal- 
culation on the non-distorted kagome lattice. The mag- 
netic intensity distribution in the scattering plane is more 
isotropic in the distorted case, consistently with the ex- 
perimental data. A second model, describing an algebraic 
spin liquid at T — 0, was derived using classical isotropic 
infinite-component spin vectors following the method of 
ref. I22T The Q distribution of magnetic intensity was cal- 
culated at a finite temperature damping the dipolar-like 
long-range spin-spin correlation by an exponentially de- 
caying term, hence smoothing the sharp features of the 
T = K calculation. The main features of the mea- 
sured Q-distribution are reproduced by both models : a 
first broad maximum localized about |Q|/|a*| = 1 from 
the origin and a second one, less well-defined, around 
|Q|/|a*| = 2.6. As expected from a static calculation, the 
best agreement between the measured and calculated Q- 
distributions is obtained at the lowest energy, 0.5 meV, at 
least concerning the position of the first maximum which 
essentially reflects the first neighbors antiferromagnetic 
correlations (cf. fig. GJ). 

The measured Q-distribution is consistent with a dis- 
order phase with specific magnetic correlations, the de- 
generacy of which would result in apparent very short 
correlation lengths. Crudely enlightened by the mod- 
els, its main features are the broad ring shape scatter- 
ing about |Q|/|a*| = 1, the scattering empty Brillouin 
zone around the origin, and the scattering minima in 
Brillouin zones around certain reciprocal lattice points 
(cf. fig. El . They can be indeed understood as ex- 
pressing the coexistence of many configurations at the 
same energy with frustration constrained antiferromag- 
netic correlations. A similar interpretation was invoked 
for the Yo.97Sco.o3Mn2 spin-liquid, the short-range mag- 
netic correlations of which would originate from 4-site 
collective spin singlets associated with the tetrahedral 
building unit of the pyrochlore lattice. Note also the 
similarity of the Q-scans recorded in both Nd-langasite 
and Yo.97Sc .o3Mn 2 .15] which can be related to the fact 
that the kagome lattice is obtained through a cut per- 
pendicular to the cube diagonal of the pyrochlore lattice. 
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A basic difference is that no dispersion of the magnetic 
signal was observed in Yo.97Sco.o3Mn2, owing to a much 
shorter lifetime of the spin correlations, « 7 10 -14 s, so 
that the modes either have no time to propagate or are 
drown within the soft modes of the ground manifold. 

Spin waves can propagate in locally ordered regions 
if the system retains a sufficient temporal and spatial 
stiffness. In a disordered medium, the spin waves life- 
time is reduced, not only through thermal fluctuations 
(i~th) HH but also through the time scale on which the 
medium fluctuates (autocorrelation time r a ) H, both 
of the order of U/ksT. The mixing of the excitations 
built on a degenerate ground state manifold further re- 

1/2 

duces their lifetime to approximately %/ [J&bT] ' (see 
ref. 0). But at the temperature we are interested in, 
we have r t h ~ T a ~ 3.8 10 -12 s and a measured life- 
time r ~ 1.4 10 -12 s. In this regime, the reduction 
due to the degeneracy of the ground state is small and 
the dynamics may still be seen as conventional e.g, the 
spin wave spectrum should be built as a superposition 
of all excitations developing over each magnetic configu- 
ration of the ground state manifold. As a consequence, 
the density of quasi static excitations should be much 
larger than the one at finite energy because of the soft 
modes. Assuming a lorentzian shape with a HWHM 
~ 0.2 meV (from r a ) for their description, it is reasonable 
to state that these quasi-static excitations also contribute 
to the scattering function at finite energy transfer. This 
would explain why the Q-dependence of the scattering 
is well accounted for by the calculated static scatter- 
ing intensities, especially at the lowest energy transfer 
(0.5 meV). Quasi static contributions from soft modes 
have a smaller weight at higher energies and propagat- 
ing modes become dominant, as evidenced from the dis- 
persion of the position of the first maximum in the Q- 
distribution. Now, a spatial propagation is expected for 
excited modes with wavelength smaller than the local 
order correlation length. The one deduced from the Q- 
dependence of the scattering (£ ~ 2 A) seems too short 
to produce the observed dispersion. This apparent con- 
tradiction can be lifted considering that the correlation 
length associated to each configuration of the ground 
state manifold would be large enough to allow propaga- 
tion of excited modes, whereas the measured one reflects 
the distribution of close configurations. The conjugate 
observation of a very broad signal associated to a dis- 
persion of its maximum could thus be a consequence of 
the high degeneracy of magnetic states in this frustrated 
compound. 

In conclusion, signatures of a cooperative paramagnet 
were evidenced in Nd3Ga5SiOi4, in particular a peculiar 
Q distribution of the magnetic correlations that would 
emerge from a highly degenerate manifold of configura- 
tions. Although the unusual dispersion of the magnetic 
correlations and the role of the single-ion magnetocrys- 
talline anisotropy have to be further clarified, the com- 



pound provides with the first measurements of the spatial 
and temporal magnetic correlations in a spin liquid phase 
of a geometrically frustrated kagome antiferromagnet. 
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